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ABSTRACT 

We report detection of continuum emission at A = 850 and 450 /im from 
disks around four Classical T Tauri stars in the MBM 12 (L1457) young associ- 
ation. Using a simple model we infer masses of 0.0014-0.012 M Q for the disk of 
LkHct 263 ABC, 0.005-0.021 M for S18 ABab, 0.03-0.18 M for LkHa 264 A, 
and 0.023-0.23 M Q for LkH« 262. The disk mass found for LkHa 263 ABC 
is consistent with the 0.0018 M inferred from the scattered light image of the 
edge-on disk around component C. Comparison to earlier 13 CO line observations 
indicates CO depletion by up to a factor 300 with respect to dark-cloud values. 
The spectral energy distributions (SED) suggest grain growth, possibly to sizes 
of a few hundred /im, but our spatially unresolved data cannot rule out opacity 
as an explanation for the SED shape. Our observations show that these T Tauri 
stars are still surrounded by significant reservoirs of cold material at an age of 
1-5 Myr. We conclude that the observed differences in disk mass are likely ex- 
plained by binary separation affecting the initial value. With available accretion 
rate estimates we find that our data are consistent with theoretical expectations 
for viscously evolving disks having decreased their masses by ~ 30%. 

Subject headings: circumstellar matter — ISM: clouds — ISM: individual (MBM 12) 
— stars: formation — stars: pre-main-sequence 
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1. Introduction 

The majority of low-mass stars emerge from their parental clouds surrounded by disks of 
0.001-0.3 M Q (Beckwith & Sargent 1996). At ages of a few Myr, these disks appear to evolve 
rapidly from optically thick at near- and mid-infrared and detectable at (sub) millimeter 
wavelengths, to undetectable at all wavelengths (Robberto et al. 2000; Jayawardhana et al. 
1999; Carpenter 2002; Haisch et al. 2001; Luhman 2001; Hartmann 2001; Duvert et al. 2000). 
The physics behind this transition and its timescale holds clues about the planet formation 
process (Clarke et al. 2001; Matsuyama et al. 2003a,b; Armitage et al. 2003). The previous 
references focus on inner disk material traced by infrared excess. Much less is known about 
colder material further from the star, even though this encompasses the bulk of the mass. 
This Letter investigates the presence of cold material around several members of the MBM 12 
young association. 

Judging from the relative occurrences of K- (~ 20%) and L- and N-band (71 ± 32%) 
infrared excess, the MBM 12 (L1457) young association is suspected to be at the very stage 
where disks start to disappear (Luhman 2001; Jayawardhana et al. 2001). After several 
unsuccessful attempts (Pound 1996, Hogerheijde et al. 2002), Itoh et al. (2003) recently 
reported detection at 1 and 2 mm of continuum emission around two or possibly three Clas- 
sical T Tauri stars in this 1-5 Myr old association (Luhman 2001), indicating the presence 
of ~ 0.05 M of cold material around each object. This Letter increases the number of 
detections of cold dust to four objects (§3), including the recently identified edge-on disk 
source LkHa 263 C (Jayawardhana et al. 2002). By extending the wavelength coverage into 
the submillimeter we can fit the spectral energy distributions (SEDs), and gain more robust 
disk-mass estimates and insight into grain growth (§4). The Letter concludes with a discus- 
sion of the inferred mass range in terms of multiplicity (Jayawardhana et al. 2002; Chauvin 
et al. 2002; Brandeker et al. 2003) and disk-dispersal models (§5). 



2. Observations 

The observations were obtained with the Submillimeter Common User Bolometer Array 
(SCUBA) (Holland et al. 1999) on the James Clerk Maxwell Telescope 5 (JCMT) on 2002 De- 
cember 12 under excellent weather conditions. Typical opacities at 225 GHz were 0.04-0.06. 
We obtained photometry of four systems: LkHa 262, the triple LkHa 263 (ABC), LkHa 264 



5 The JCMT is operated by the Joint Astronomy Centre in Hilo, Hawaii, on behalf of the parent orga- 
nizations PPARC in the United Kingdom, the National Research Council of Canada and The Netherlands 
Organization for Scientific Research. 
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(A), and the triple S18 (ABab); Table 1 lists coordinates and observing details. The em- 
ployed two-bolometer photometry technique allowed for increased observing efficiency, with a 
chop throw of 52'.'4 in Naysmith coordinates. Individual integrations of 30 s were repeated for 
the totals listed in the table. To confirm the unresolved nature of the emission we obtained 
a 64-point jiggle map centered on LkHa 262 also containing LkHa 263 (ABC). While pro- 
viding spatial information, the jiggle map is less sensitive than the single-point photometry 
data. 

The standards HL Tau, CRL 618, and Uranus provided focus checks and flux calibra- 
tions; the nearby quasar 0235+164 served as pointing source every ~ 2 hours. Pointing 
accuracy was good with excursions of less than a few arcsec. However, uncorrected pointing 
errors can still affect the photometry in the 8" beam at 450 /xm. In spite of a 3" dither 
included in the photometry, 450 /im fluxes of LkHa 262 from photometry are lower by 60% 
than those from the jiggle map (§3). We conclude that uncorrected pointing offsets and cal- 
ibration uncertainties at short wavelengths due to the imperfect beam shape of the JCMT 
resulted in 450 /im photometry results that are strict lower limits to the actual source flux. 
We include a +60% error in the uncertainty of the reported 450 /im results. The 850 /im 
photometry and the jiggle maps are unaffected. 



3. Results 

All four objects, containing a total of eight (known) stars, show emission at 850 and 
450 /iiii (Table 2). The jiggle map (Fig. 1) indicates that the emission is unresolved and 
confined to the source position, and is not structure in the cloud. Because of the higher 
noise level, only LkHa 262 is detected in Fig. 1 while LkHa 263 (ABC) remain undetected. 
The separation of 17" between LkHa 262 and 263 (ABC) is large enough that < 1-2 mJy 
spill over at both wavelengths is expected, based on archival beam profiles. 

High signal-to-noise photometry and jiggle-map data of LkHa 262 are consistent at 
850 /im but discrepant at 450 /im with respective fluxes of 166.2 mJy and 263.1 mJy. Un- 
corrected pointing and calibration errors at 450 /im are likely to blame (§2); the jiggle-map 
flux is extracted after Gaussian profile fitting and does not suffer from pointing offsets. The 
reconcile the measurements, we include a +60% error bar in Table 2 and stress that the 
450 /im photometry values are strict lower limits. The 850-450 /im spectral indices a conse- 
quently contain a large uncertainty (Table 2). At their high end, corresponding to the high 
end of allowed 450 /im fluxes, the indices (a =1.5-2.5) are consistent with emission from 
cool (< 30 K) and coagulated dust grains. The spectral index a = (3 + 7, where the index 
of the dust emissivity (5 decreases when grains coagulate (e.g., Ossenkopf & Henning 1994; 
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Pollack et al. 1994) and where the slope of the Planck function 7 falls below the value of 2.0 
outside the Rayleigh- Jeans limit when T < hu/k = 32 K at 450 fxm. 

4. Disk models and masses 

Fig. 2 plots our data and values from Pound (1996), Hogerheijde et al. (2002), and Itoh 
et al. (2003); the +60% error bars to the 450 /im photometry are explicitly included. From 
these SEDs we can infer the mass of cold material. The 850 fxm fluxes trace the absolute 
amount of material, while the SEDs help to constrain important model parameters. We 
choose the flared-disk model of Chiang & Goldreich (1997) to describe the disks. We use the 
disk mass as our only free parameter. The disk temperature distribution is fixed by scaling 
to the stellar luminosity oc (L^/Lo) -0 ' 25 , where Lq = 1.42 L Q (Chiang & Goldreich 1997). 
Source luminosities L± are from Luhman (2001). We neglect any changes in disk structure 
due to changes in mass and temperature. Following recent reassessments of the distance 
to MBM 12 (Luhman 2001; Andersson et al. 2002) we adopt 275 ± 60 pc. We assume an 
average inclination of 60°. Because most of the flux comes from optically thin regions, only 
near edge-on orientations would change our results significantly. 

The final model parameter is the dust emissivity k v : its absolute value and its variation 
with wavelength. Often, k v is parameterized as k v = kq{uqY with k = 10. cm 2 g" 1 (dust) 
and u = 10 12 Hz (Hildebrand 1983). From disk SED fitting, low values of (3 are commonly 
found, 0-1 (Beckwith & Sargent 1991) and interpreted as evidence for grain growth; Beck- 
with et al. (2000) point out that grain mineralogy can affect (5 and explain the low values. 
More accurate descriptions of k v (v) follow from calculations of grain growth in astrophysi- 
cal environments employing realistic mineralogies (e.g., Ossenkopf & Henning 1994; Pollack 
et al. 1994; Beckwith et al. 2000). These models circumvent having to choose a value for k , 
but still show show significant variation. We adopt a dual approach: we fit the SEDs with 
a parameterized k v to constrain f3 and learn about possible grain growth, and then fit the 
SEDs with several physical grain-growth models to derive accurate disk masses with realistic 
uncertainties. 

With k v = Ko{v Q y, x 2 minimization of the SED fit yields (5 = 0.5 ± 0.3 for LkHa 262 
and 0.8 ± 0.5 for LkHa 264; the error bars for LkHa 263 and S18 are too large to constrain 
p. For this fit we have used all detections and upper limits at wavelengths > 450 /im. At 
the wavelengths of the infrared N-band and shorter, the flux is dominated by very small 
amounts of hot material and may contain contributions from silicate emission; both depend 
strongly on the model parameters. This does not affect our results: < 1% of the 450 /im 
flux originates from the < 4 AU region that contributes > 96% of the 10 /im emission. The 
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inferred values (3 — 0.5 ± 0.3 and 0.8 ± 0.5 are similar to those commonly found in T Tauri 
disks, and we interpret them as indications for grain growth in the MBM 12 disks. 

Few of the grain-growth models of Ossenkopf & Henning (1994) and Pollack et al. 
(1994) have f3 < 1. Of the models by Ossenkopf & Henning, only strongly coagulated grains 
without ice mantles show such low (5. The authors note that such grains are unlikely in disk 
environments and we follow their advice and avoid this particular class. Of the models by 
Pollack et al., f3 < 1 is only found for average grain sizes > 1 cm. Table 2 lists the range of 
disk masses found by x 2 minimization. These vary from 0.0014-0.012 M Q for LkHa 263 to 
0.023-0.23 M Q for LkHa 262. The allowed range for each source reflects the uncertainty in 
the dust opacity models. These these ranges overlap, but differences in disk mass between 
sources are robust if we assume similar dust properties in all objects. 

Our modeling makes no assumptions about the opacity. We find that < 15% of the 
450 /iin flux and < 7% at 1.3 mm originates from opaque regions (r > 3). Only much 
smaller disks (< 30 AU) have ~ 85% of the 450 /iin flux coming from optically thick regions 
(45% at 1.3 mm). For such disks the SED slope can be explained by opacity instead of grain 
growth, and submillimeter fluxes no longer trace mass. Only spatially resolved data can 
settle this issue. 

Jayawardhana et al. (2002) infer 0.0018 M for LkHa 263 C from modeling the scat- 
tered light, depending on the assumed dust properties. Our values of 0.0014-0.012 M are 
consistent with their findings. An unknown fraction of our inferred mass may reside in disks 
around the companions A and B. Spatially resolved data and detailed modeling are required 
to further characterize the disk of LkHa 263 C. 

13 CO 2-1 measurements (Hogerheijde et al. 2002) indicate that CO is depleted with 
respect to the dark cloud value of CO/H 2 =2 x 10" 4 by factors up to 300. Large depletion 
levels are found in disks around T Tauri stars (e.g., Dutrey et al. 1994, 1996) and expected 
theoretically (e.g., Aikawa & Herbst 1999; Willacy & Langer 2000). 

5. Nature or nurture 

Although our sample is small, with only eight stars in four unresolved systems, interest- 
ing conclusions can be reached about the disk evolution. Our detections show that at least 
the Classical T Tauri stars in MBM 12 (1-5 Myr) have disks in the same mass range as the 
younger (~ 1 Myr) Taurus and p Ophiuchus regions (0.001-0.3 M ; Beckwith & Sargent 
1996). While Carpenter (2002) shows that 0/95 stars in IC 348 (~ 2 Myr) have disk masses 
> 0.025 M , two out of four of our MBM 12 systems do (or > 17% of all twelve MBM 12 
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systems); in Taurus, a comparable fraction of 14% of stars have disks > 0.025 M Q . Luhman 
(2001) argues that the K- and L-band excess suggest significant disk dispersal in MBM 12. 
Our sample does not directly address this issue because 3/4 objects have K-band excess (S18 
does not). 

Jensen et al. (1996) showed that binaries with separations < 50-100 AU have reduced 
disk mass or no disks, compared to wider binaries and single stars. Our data follow that 
trend, taking into account that the relative disk masses are robust (§4): LkHa 262: 15'/3 = 
4208 AU and 0.023-0.23 M ; LkHa 264 A: 9"2 = 2530 AU and 0.03-0.18 M Q ; S18 ABab: 
0"75 = 206 AU (A-B), 0"063 = 17 AU (Ba-Bb) and 0.005-0.025 M ; and LkHa 263 ABC: 
4'/l = 1128 AU (AB-C), 0"41 = 113 AU (A-B) and 0.0014-0.012 M (binary separations 
from Jayawardhana et al. 2002, Chauvin et al. 2002 and Brandeker et al. 2003). This 
suggests that the presence of a close binary (environment) determines the total amount of 
disk material. 

If environment determines (initial) disk mass, can we still infer anything about disk 
dispersal from the amount of material? In a viscouly evolving disk, disk mass and the 
accretion rate onto the central star are linked (Clarke et al. 2001; Matsuyama et al. 2003a). 
Fig. 3 plots 'evolutionary tracks' of disks with initial masses of 0.001, 0.01, and 0.1 M Q 
around a 0.5 M star, following Matsuyama et al. (2003a). This model assumes a standard 
a-disk description with a = 10~ 3 , a central star with = 0.5 M and /?* = 2 R , and a 
constant ionizing flux from the star with time- dependent contribution from accretion. The 
disk removal rate only depends weakly on the stellar parameters, and any uncertainties are 
insignificant compared to our observational error bars. 

Muzerolle et al. (in prep.) report estimates of mass accretion rates onto three of our 
objects from Br7 line measurements, using the calibration of Muzerolle et al. (1998). They 
find upper limits of (4-5) xlO -9 M yr _1 for LkHa 263 and 264, and a rate of (2 ± 1) x 
10~ 9 M yr _1 for S18. These rates are comparable to those inferred for IC 348 (Liu et al. 
2003) and at the low end of the range found in Taurus, which has a median of 10~ 8 M yr -1 . 
The errors include uncertainties in the observations and the stellar mass, assumed to be 
0.5 M , but not the distance. Because both accretion rate and disk mass depend on the 
distance squared, its uncertainty does not enter into the values of Fig. 3; M acc also depends 
on M+/R+, which we assume constant. 

The disk mass and accretion rate of S18 are consistent with the dispersal model for 
the age range of MBM 12 (1-5 Myr; shaded in Fig. 3), suggesting that the object has lost 
one-third of its initial disk mass of ~ 0.016 M . The upper limit on the accretion rate 
of LkHa 264 also agrees with the model at 1-5 Myr. The agreement improves if the disk 
mass is decreased somewhat, which is expected when dust grains have indeed grown to a 
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few hundred jum as suggested by the SED. Grain sizes comparable to observing wavelength 
are efficient emitters, reducing the required mass to fit the observed flux. The limits on 
LkHa 263 do not provide useful constraints. 

In summary, we conclude that the Classical T Tauri stars in MBM 12 still have significant 
reservoirs of cold dust in circumstellar disk. There are indications for grain growth up to 
several hundred /im in these disks, but spatially resolved observations are required to rule 
out opacity as an explanation for the flat spectral slopes. And while differences in disk 
mass are likely dominated by environment (binary separation; 'nature'), available accretion 
rates and detected disk masses are consistent with a disk dispersal scenario ('nurture') with 
one-third of the mass already lost. 

We thank our TSSs J. Kemp and J. Hoge for excellent support during our observations. 
The staff of the JCMT and the JAC are thanked for their hospitality. RJ acknowledges 
support in part through NASA Origins Grant NAG5-11905. I.M. acknowledges support 
from the National Reserach Council, Canada. The referee is thanked for a careful reading 
of our manuscript and insightful comments. 



- 8- 



REFERENCES 

Aikawa, Y. & Herbst, E. 1999, A&A, 351, 233 

Andersson, B.-G., Idzi, R., Uomoto, A., Wannier, P. G., Chen, B., & Jorgensen, A. M. 2002, 
AJ, 124, 2164 

Armitage, P. J., Clarke, C. J., & Palla, F. 2003, MNRAS, in press 

Beckwith, S. V. W., Henning, T., & Nakagawa, Y. 2000, Protostars and Planets IV, 533 
Beckwith, S. V. W. & Sargent, A. I. 1991, ApJ, 381, 250 
-. 1996, Nature, 383, 139 

Brandeker, A., Jayawardhana, R., & Najita, J. 2003, AJ, submitted 
Carpenter, J. M. 2002, AJ, 124, 1593 

Chauvin, G., Menard, F., Fusco, T., Lagrange, A.-M., Beuzit, J.-L., Mouillet, D., & 
Augereau, J.-C. 2002, A&A, 394, 949 

Chiang, E. I. & Goldreich, P. 1997, ApJ, 490, 368 

Clarke, C. J., Gendrin, A., & Sotomayor, M. 2001, MNRAS, 328, 485 

Dutrey, A., Guilloteau, S., Duvert, G., Prato, L., Simon, M., Schuster, K., & Menard, F. 
1996, A&A, 309, 493 

Dutrey, A., Guilloteau, S., & Simon, M. 1994, A&A, 286, 149 

Duvert, G., Guilloteau, S., Menard, F., Simon, M., & Dutrey, A. 2000, A&A, 355, 165 
Haisch, K. E., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, L153 
Hartmann, L. 2001, AJ, 121, 1030 
Hildebrand, R. H. 1983, QJRAS, 24, 267 

Hogerheijde, M. R., Jayawardhana, R., Johnstone, D., Blake, G. A., & Kessler, J. E. 2002, 
AJ, 124, 3387 

Holland, W. S., Robson, E. I., Gear, W. K., Cunningham, C. R., Lightfoot, J. F., Jenness, 
T., Ivison, R. J., Stevens, J. A., Ade, P. A. R., Griffin, M. J., Duncan, W. D., Murphy, 
J. A., & Naylor, D. A. 1999, MNRAS, 303, 659 



- 9- 



Itoh, Y., Sugitani, K., Fukuda, N., Nakanishi, K., Ogura, K., Tamura, M., Marui, K., Fujita, 
K., Oasa, Y., & Fukagawa, M. 2003, ApJ, 586, L141 

Jayawardhana, R., Hartmann, L., Fazio, G., Fisher, R. S., Telesco, C. M., & Pifia, R. K. 
1999, ApJ, 521, L129 

Jayawardhana, R., Luhman, K. L., D'Alessio, P., & Stauffer, J. R. 2002, ApJ, 571, L51 

Jayawardhana, R., Wolk, S. J., Barrado y Navascues, D., Telesco, C. M., & Hearty, T. J. 
2001, ApJ, 550, L197 

Jensen, E. L. N., Mathieu, R. D., & Fuller, G. A. 1996, ApJ, 458, 312 

Liu, M. C, Najita, J., & Tokunaga, A. T. 2003, ApJ, 585, 372 

Luhman, K. L. 2001, ApJ, 560, 287 

Matsuyama, I., Johnstone, D., & Hartmann, L. 2003a, ApJ, 582, 893 
Matsuyama, L, Johnstone, D., & Murray, N. 2003b, ApJ, 585, L143 
Muzerolle, J., Hartmann, L., & Calvet, N. 1998, AJ, 116, 2965 
Ossenkopf, V. & Henning, T. 1994, A&A, 291, 943 

Pollack, J. B., Hollenbach, D., Beckwith, S. V. W., Simonelli, D. P., Roush, T., & Fong, W. 
1994, ApJ, 421, 615 

Pound, M. W. 1996, ApJ, 457, L35 

Robberto, M., Beckwith, S. V. W., & Herbst, T. M. 2000, in Star Formation 1999, ed. 
T. Nakamoto (Nobeyama: Nobeyama Radio Observatory), 231 

Willacy, K. & Langer, W. D. 2000, ApJ, 544, 903 



This preprint was prepared with the AAS IATgX macros v5.0. 



-10- 



Fig. 1. — Jiggle maps at (a) 850 jum and (b) 450 /im containing four (known) young stars 
indicated by symbols: LkHa 263 C, LkHa 263 A+B, and LkHa 262 (left to right). The 
emission around LkHa 262 is unresolved. Emission from LkHa 263 ABC is undetected at 
the noise level of 7 mJy beam -1 (850 /im) and 22 mJy beam -1 (450 /im). 

Fig. 2. — Spectral energy distributions (SEDs) of the four sources detected with SCUBA. 
Filled symbols are SCUBA data, open symbols data from Itoh et al. (2003), star symbols 
show N-band excess (Jayawardhana et al. 2002), and 2a upper limits are from Pound (1996), 
Hogerheijde et al. (2002), and Itoh et al. (2003). The broad shaded curves show the model 
SEDs and their uncertainty due to variations in adopted dust emissivities. Emission at 
< 20 /im originates from small amount of warm dust and is highly model dependent. 

Fig. 3. — 'Evolutionary tracks' of viscous accretion disks from Matsuyama et al. (2003a) 
(solid curves), for a stellar mass of 0.5 M and initial disk masses of 0.1, 0.01, and 0.001 
M . Dashed lines are isochrones at 1 x 10 4 , 4 x 10 4 , 1.6 x 10 5 , ... yr. At these ages, the disks 
are reduced to 99%, 97%, 91%, 75%, 50%, and 27% of their initial mass. The shaded area 
is the inferred age range of 1-5 Myr of MBM 12. Upper limits are 2a values. 
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Table 1. Observations 





a(2000) 


5(2000) 


Observing 


Integration 


Source 


(h m s) 


(° ' ") 


Mode 


Time (s) 


LkHa262 


02 56 07.9 


20 03 25 


photom 


1200 








jiggle 


4800 


LkHa263 ABC a 


02 56 08.7 


20 03 41 


photom 


4800 


LkHa264 A 


02 56 37.5 


20 05 38 


photom 


1200 


S18 ABab 


03 02 21.1 


17 10 35 


photom 


1200 



a Centered on LkHa263 C, 4"1 from AB. 



Table 2. Results 





F(850 fxm) 


F(450 iim) 


Spectral Disk Mass 


Source 


(mJy) 


(mJy) 


Index a (1(T 2 M Q ) 



LkHa 262 (photom) 97.2 ± 3.4 166.21?^ 0.67-1.64 0.023-0.23 a 

LkHe* 262 (jiggle) 95.4 ± 7 263.1 ±22 1.35-1.84 

LkHa 263ABC 7.5 ± 1.8 26.7^ 6 2 ° 1.16-3.17 0.0014-0.012 a 

LkHct 264A 136.7 ± 3.4 214.1±^ s 9 5 0.58-1.48 0.03-0. 18 a 

S18ABab 21.6 ±3.3 65.2+^ 1.20-2.74 0.005-0.025 a 



a The lower mass limits corresponds to dust opacities from Ossenkopf & Henning (1994) 
with ice mantles and > 10 6 yr of grain growth; the upper mass limits to opacities from 
Pollack et al. (1994) with grain sizes ~ 1 cm. 
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